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Abstract 

We investigate the effects of low energy heavy nuclei from the supernova 
explosion on nearby interstellar space. In addition to the ionization and heating 
of the Gum Nebula, these particles may produce detectable fluxes of X-rays and 
gamma rays, both as continuum radiation and line emission. 


Introduction 

The possibility that the ionization of the Gum Nebula may have been produced 
by energetic charged particles from the supernova Vela X has been suggested by 
Ramaty et al. (1971). The principal features of this model are the existence of a 
prolonged source of ionization and heating and the possibility of direct detection 
of hard photon emissions from contemporary fluxes of energetic particles in the 
nebula. In the present paper, we shall summarize the model and discuss its main 
features. We shall also present an estimate of gamma-ray line emission from the 
interaction of energetic nuclei with the ambient gas and we shall discuss some 
additional ideas regarding the maintenance of a possible high temperature in the 
filaments of the nebula. 


Particle Propagation and Energy Deposition 

A necessary requirement of a cosmic-ray model of the Gum Nebula is that 
in a time less than or equal to the age of the supernova remnant the particles 
should be capable of both reaching the outer edges of the nebula and depositing 
a significant fraction of their energy as ionization loss to the ambient gas. 

A lower limit on the velocity c (3 of the ionizing nuclei can be obtained by 
assuming rectilinear propagation of particles in the nebula, i.e. : 

J3 > L/ct (1) 
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where L and t are the nebular radius and age respectively. If we adopt an age of 
about 10 4 years, as deduced from the period and rate of change of period of PSR 
0833-45 (Reichley, Downs and Morris 1970), and a nebular radius of about 360 pc, 
as obtained (Alexander et al. 1971) from La absorption data and a local hydrogen 
density of 0.24 cm -3 , the lower limit on /3 is about 0.1 corresponding to an energy 
per nucleon e ^ 5 MeV/nucleon. It should be remembered that this lower limit 
is subject to considerable uncertainties in both t and L. In particular, if the age 
is larger than 10 4 years (Shklovsky 1970), the Gum Nebula may have been pro- 
duced by even lower energy particles . 

Cosmic-ray particles in general do not propagate along rectilinear paths but 
follow magnetic field lines and are scattered by irregularities in the field. The 
mean distance traversed in such diffusive motion may be approximated by 

L 2: (2\c/?t) 1/2 (2) 

where A. is the diffusion mean free path. For L - 360 pc and t - 10 4 years, 

A/3 -20 pc. If /3 = 0.15 (e = 10 MeV/nucleon), A - 130 pc. This value is quite 
large but not inconsistent with recent ideas on cosmic ray propagation (Ramaty 
and Lingenfelter 1971). In particular, if cosmic rays propagate by compound 
diffusion (Lingenfelter, Ramaty and Fisk 1971) or if cosmic rays are confined 
to the galactic disk by reflecting boundaries, both the low anisotropy and abun- 
dances of fragmentation products of galactic cosmic rays are not inconsistent 
with the long diffusion mean free path required for the ionization of the Gum 
Nebula. Finally, a relatively scatter-free motion in the nebula is also supported 
by the low values of both the ambient gas density and magnetic field in this 
region of interstellar space. 

The ionization loss rate and the time for the deposition of half the energy of 
a nonrelativistic particle of mass number A and nuclear charge Ze in neutral 
hydrogen of density n H are given by 


— * 1.46 x 10' 12 Z 2 ,/A n„ e ' 0 - 3 MeV/nucl. sec " 1 
dt eff H 


(3) 


and 


t l/2 ~ 104 ’ e1 ' 3 A^eff n H> y ears 


( 4 ) 


where e and n„ are in MeV/nucleon and cm -3 , respectively. As a result of 

ri 

electron capture during the slowing down process the effective charge Z eff 
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becomes smaller than the nuclear charge Z. Z /Z may be approximated by 
(Pierce and Blann 1968) 


Z eff /Z ^ 1 - exp (- 130 yfl/Z 2/3 > 


( 5 ) 


Since even for Fe 56 at 10 MeV/nucleon Z f /z - 0.9, we shall use Z instead of 
Z ef f in the evaluation of equations (3) and (4). 

In Table 1, n H t 1/2 is given for various 10 MeV/nucleon particles. As can 
be seen, in 10 4 years these particles will deposit a major fraction of their 
energy only if n H » 1 cm -3 . The average electron density <n e > in the Gum 
Nebula is only 0.16 cm -3 , but the density in filaments is much higher, n e - 
X <n g > where X ^65 is the clumping factor in the nebula (Brandt et al. , 1971). 

A possible model, therefore, is one in which initially the nebular region con- 
sisted of neutral clouds or clumps with hydrogen densities approximately equal 
to the presently observed electron densities in filaments. The particles ejected 
from the supernova propagate in a scatter- free manner in the intercloud medium 
and subsequently penetrate clouds and lose their energy in a time scale shorter 
than the trapping time in the cloud. A possible trapping mechanism, based on 
wave-particle interactions, has been discussed by Ramaty et al. (1971). 

If the initial hydrogen density in clouds is equal to the present electron den- 
sity of about 10 cm -3 , it is very unlikely that low energy protons have ionized 
the Gum Nebula, since, as can be seen from Table 1, it would take about 2 x 10 4 
years for these particles to lose half their energy. On the other hand, if the 
ionization were produced by heavy nuclei such as Si and Fe, these particles 
would deposit sizable fractions of their energy in clouds on a time scale of 2000 
to 3000 years leaving about 7000 to 8000 years for scatter-free propagation in 
the intercloud medium. This is consistent with the estimates of the transit 
times to the edges of the nebula given above. 


Table 1 



n H t i/2 < cm ~ 3 years) 

Hi 

200000 

C 12 

67000 

Si 28 

29000 

Fe 56 

16000 
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A specific consequence of this model is that the degree of ionization of the 
intercloud medium is much smaller than that of the filaments. This is not the 
case for models in which the ionization is produced by photons since these lead 
to the complete ionization of both the clouds and intercloud medium out to an 
appropriate Stromgren radius. 


Mass and Energy Requirements of Supernova Ejecta 

A fast particle of energy £ 0.1 MeV/nucleon expends some 36 eV to produce 
an ion pair (Bethe and Ashkin 1953). Therefore, if the nebula contains N e ~ 

2 x 10 62 electrons, as suggested by Brandt et al. (1971), the particles ejected 
by the supernova must lose at least 1.2 x 10 52 ergs. For the present value of 
N e , this is a lower limit, since a fraction of the energetic particles may slow 
down in an ionized medium, thereby converting their energy directly to heat, 
and the total number of electrons produced could be larger than N e since some 
recombination in the filaments may have already taken place. 

The velocity distribution of the ejected material from a supernova explosion 
may be approximated by (Colgate and McKee 1969) 


iv. (1 - F) 1/3 : 


0.75 £ F £ F n = 3/7 


V- - 


( 6 ) 


0.67 F” 1/4 ; F < F 0 


where F is the external mass fraction. The corresponding number spectrum of 
nuclei of mass number A in the nonrelativistic region is given by 


r, 


6 M ei ( A ) 1 

N (e) = - -1 1 _ J 

W 7 Am e ^ 


O/e 0 ) 


3/2 . 


e - 


(V e ) 2 ; e> 


(?) 


Upon integration, the total kinetic energy W of the ejected material is 


M . 

*4— «c 

5 m L 


( 8 ) 


where m p is the proton mass. 
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Since W > 1.2 x 10 52 ergs, the ejected mass in units of a solar mass and 
the characteristic ejection energy in MeV/nucleon must obey 

M . 

— — - « 0 £ 5 MeV/nucleon (9) 


As discussed above, « 0 cannot be smaller than about 5 MeV/nucleon, since 
otherwise the particles will not reach the edge of the nebula in the available time 
of about 10 4 years. On the other hand, e 0 cannot be much larger than about 20 
MeV/nucleon, since, as can be seen from Table 1 and equation (4), the ejected 
nuclei, even if they are all Fe 56 , will not have enough time to deposit their 
energy in clouds. For the range 5 < e. Q < 20 MeV/nucleon, the lower limits on 
the ejected mass range from 0.25 M Q to 1 M 0 . For the efficient ionization 
of the nebula, the ejected material must consist of nuclei heavier than at least 
C 12 . 


We now proceed to examine the observational consequences of the model. 


Hard Photon Emissions by Energetic Particles in the Nebula 

X-ray emission by bremsstrahlung is a necessary consequence of the 
interaction of fast nuclei with interstellar gas (Boldt and Serlemitsos 1969). 

The photon energy spectrum that would result from the particle distribution 
given in equation (7) is essentially flat with an end point hv 0 at (m e /m p )e 0 , 
where m e is the electron mass. 

Let M/(Am ) nuclei of mass number A and atomic number Z be currently 
trapped in clouds of density n e . The average bremsstrahlung volume emissivity 
in the Gum Nebula (independent of the ionization state of the gas) is given by 


7)--— Z 2 a Q m c 2 — <c/3> (10) 

77 0 6 V Am V V ’ 

p 

where cr Q = 6.7 x 10" 25 cm 2 , a = 1/137, V is the volume of the Gum Nebula, and 
<c/5> is the particle velocity averaged over the energy spectrum of the par- 
ticles. For 0} = 0.15, Z 2 /A = 8.7 (average of Si 28 and Fe 56 in a ratio of 2 to 
1), n e = 10 cm~ 3 and V = 1.35 x 10 63 cm 3 (Brandt et al. , 1971) we get 


77 — ■ 9 x 10" 28 


M 

M e 


erg cm 3 sec* 1 


( 11 ) 
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According to the geometry of the Gum Nebula described by Brandt et al. 
(1971), we expect it to appear as a disk source of X-rays, extended in galactic 
longitude (A/ - 90°) and relatively thin in galactic latitude (Ab - 7°). Collecting 
radiation from an interval of length L penetrating the nebula, the intensity 
viewed by a detector of aperture (Ab) 0 < Ab would be expressed by 

dl/d't = (Ab) 0 /(47r) t?L erg cm" 2 sec -1 rad -1 (12) 


Using a detector of (Ab) 0 = 4°, viewing at l - (260 ± 40)°, which spans the 
Gum Nebula, Cooke, Griffiths and Pounds (1969) observed an enhanced X-ray 
emission associated with the galactic disk of about 3 x 10 " 9 erg (cm 2 sec rad) -1 
above 1.4 keV, with an essentially flat energy spectrum extending to about 10 
keV. The upper limits to this flux at energies above 12.5 keV, as set by Hudson, 
Peterson, and Schwartz (1971), provide evidence that the spectrum is character- 
ized by a break in the vicinity of about 10 keV. Using the volume emissivity as 
obtained from equation (11), L = 400 pc and (Ab) 0 = 4°, we get dl/d't = 6 x 10 " 9 
M/M 0 erg sec -1 cm" 2 rad" 1 . Since a significant part of the apparently diffuse 
background observed from the direction of the Gum Nebula may be coming from 
undetected discrete sources (Cooke and Pounds 1971), the flux of 3 x 10 " 9 erg 
cm" 2 sec" 1 rad' 1 must be considered as an upper limit. This then sets an 
upper limit M/M q < 0.5 on the mass of energetic charged particles currently 
trapped in clouds. 

In addition to bremsstrahlung X-rays, contemporary fluxes of energetic 
heavy nuclei in dense clouds may also produce detectable gamma- ray line 
emission in the few MeV range. As an example, we consider the 1.78 MeV line 
of Si 28 . The cross section cr^ for the reaction Si 28 (p, p’) si 28 * 1 78 MeV has a 
threshold at about 5 MeV, peaks at about 9 MeV and has a value of about 300 mb 
at 10 MeV (McGowan et al. 1969, 1970). 

As above we assume that M/(Am ) particles are currently trapped in 
clouds. The average volume emissivity of the Gum Nebula in line emission at 
a photon energy e. is then given by 


1 («j) =cr j ( n p /V) (2M/3Am p ) <c /3> 


(13) 


where n p equal to n e is the proton density in clouds, and the factor of 2/3 is the 
assumed fraction of silicon nuclei in the supernova ejecta. For the values of 
n = n , V and </ 3 > given above 

P ® 
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(14) 


7] (1.78 MeV) - 8 x 10" 28 — erg cm 3 sec -1 

M o 

Since present measurements of gamma rays at these energies have been made 
with omnidirectional detectors only (Vette et al. 1970), we shall make the 
simplifying assumption that the flux at earth is produced by a point source at 
the center of the Gum Nebula. The flux seen by an omnidirectional detector 
will then be given by 


F (e. ) rj (e. ) V/(4 77 r 2 e. ) (15) 

For r - 460 pc, e. = 1.78 MeV and M/M 0 < 0.5, F (1.78 MeV) < 8 x 10" 3 photon 
cm' 2 sec _1 . This flux should be compared with the omnidirectional counting 
rate of 1.2 x 10" 1 photon cm -2 sec -1 in the 1 to 2 MeV range (Vette et al. 1970). 
Since the energy of the photon- emitting nuclei is about 10 MeV/nucleon or less, 
line broadening should be less than 15%. Therefore, detectable gamma- ray line 
emission could be observed with detectors of moderate energy resolution. 

It is of some interest to speculate about line emission at 4.43 MeV from 
C 12 and 6.14 MeV fromO 16 . The cross sections for the excitation of these 
levels were plotted as a function of energy by Lingenfelter and Ramaty (1967). 

At 10 MeV/nucleon cr - 200 mb for C 12 and a - 100 mb forO 16 . If M 0 and 

ry fy W 

M 0 are the present masses of energetic carbon and oxygen nuclei in clouds the 
expected photon fluxes from the Gum Nebula are F(4.43 MeV) - 4 x 10" 2 M c /M 0 
photon cm _2 sec _1 and F(6.14 MeV) - 1.5 x 10" 2 M c /M 0 photon cm" 2 sec _1 . When 
compared with the observed counting rate of about 5.6 x 10" 2 photon cm -2 sec -1 
in the 4 to 6 MeV band (Vette et al . 1970), we see that line emission from car- 
bon and oxygen could be detected and, moreover, it could account for at least 
part of the excess radiation above the intergalactic background as calculated 
by Brecher and Morrison (1969). 

In addition to the gamma ray region, line emission should also be produced 
in the X-ray region as a result of charge exchange between energetic nuclei and 
H atoms resulting in the capture of electrons to excited states. Cascades down 
to the ground state occur, which produce the analogue of Ly a and Ly/3 emission 
for the case of capture of K-electrons. Silk and Steigman (1969) have given a 
discussion of this process, with application to low energy galactic cosmic rays. 
The contemporary flux of X-ray line emission, however, is difficult to determine, 
since unlike bremsstrahlung and nuclear collision which do not depend on the 
ionization state of the ambient matter, charge exchange takes place in neutral 
hydrogen only. Estimates of the time integrated fluxes in the 7 keV line of iron 
and 2 keV line of silicon were given by Ramaty et al. (1971). 
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Cooling of the Nebula 


Alexander et aL (1971) suggested that the temperature of the nebula as a 
whole is quite high, probably in excess of 5 x 10 4 °K. While it is possible that 
the dense filaments could at present be cooler than 10 4 °K, it is of some 
interest to examine the effects of a contemporary source of ionization and heat- 
ing on the thermal state of the filaments . 

According to Cox and Tucker (1969), for a gas of cosmic abundances in 
ionization equilibrium, the radiative loss rate at 5 x 10 4 °K is P - 2 x 10 22 
n 2 erg cm -3 sec -1 . The corresponding cooling time is t c -(3/2)n e kT/P - 
1700/n e years. For n e = 10 cm -3 , t c -170 years, much shorter than the age 
of the nebula. 

The rapid cooling of the filaments can be partially compensated by a con- 
temporary source of heat such as a large flux of energetic particles presently 
trapped in clouds. The rate of energy loss of a particle of mass number A and 
nuclear charge Ze in a plasma of density n e is approximately given by equation 
(3) with n H replaced by 4n e (Hayakawa and Kitao 1956). If M/(Amp) nuclei are 
currently trapped in clouds, the contemporary heat input is given by 


Q - 5-6 x 10“ 22 (Z 2 / A) e “ 0, 3 n — erg cm 3 sec " 1 (16) 

6 M o 

For Z 2 /A = 8.7, e = 10 MeV/nucleon, n =10 cm -3 and M/lVL < 0.5, Q < 0.6P. 
Thus, even though a contemporary flux of energetic particles, equal to the upper 
limit determined from the calculated bremsstrahlung emission and the observed 
disk component of diffuse x-rays could make a contribution to the maintenance 
of the temperature in the filaments, it is not quite sufficient to balance the radi- 
ative energy loss of the plasma. The characteristic time scale for contemporary 
heating, t h , is given by equation (4) with n H replaced by 4n . For Z 2 /A = 8.7, 

« = 10 MeV/nucleon and n e = 10 cm -3 , t 800 years. 

The radiation loss P, however, could be smaller than the value calculated 
by Cox and Tucker (1969) if the abundances of carbon and oxygen in the Gum 
Nebula are below universal abundances or if the ions C + 1 to C +3 andO +1 to 
0 +s are significantly depleted from their values corresponding to ionization 
equilibrium by the initial ionization process. If the maintenance of a high 
temperature in the filaments is achieved by the latter mechanism, the charac- 
teristic time of interest is the recombination time ofO +6 toO +s and C +4 to C +3 . 
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Up to temperatures of a few times 10 s °K dielectronic recombination for 
0 +6 and C +4 is negligible. Using Seaton’s (1959) hydrogenic approximation, the 
radiative recombination coefficient to all levels of an ion X +m is given by 
( e.g. : Allen and Dupree 1969) 


a - 2 x 10" 11 (m + l) 2 T" 1 - /2 (0.43 + — in A. + — ] cm 3 sec -1 (1^) 

V 2 \ i/3 J 

where K = 157890 (m H ) 2 /T. For n e = 10 and T = 2 x 10 s °K, the recombination 
time for bothO +6 - 0 +s and C +4 - C +3 is about 900 years. By combining this 
time scale with the contemporary heating time t h , we see that the onset of rapid 
cooling could be delayed by as much as 1500 to 2000 years. A high temperature 
in the filaments could then be consistent with a delayed mechanism for the ioni- 
zation and heating of the Gum Nebula. 


Ionization and Heating of the Interstellar Medium 

Several recent theoretical descriptions (Pikel'ner 1967; Bala Subrahmanyan 
et al. 1968; Spitzer and Tomasko 1968; Spitzer and Scott 1969; Goldsmith, Habing, 
and Field 1969) of ionization and thermal equilibria for interstellar H I regions, 
based upon current models for the interstellar gas, indicate that the rate of ion- 
ization per hydrogen atom is £ - 10" 15 (sec H atom) - J . Hjellming, Gordon, and 
Gordon (1969) find that observed pulsar dispersion measures may be best fitted 
with such a model for i = (2.5 ± 0.5) x 10" 15 (sec H atom) -1 . For an equili- 
brium situation, the rate of electron-ion recombination is the most direct meas- 
ure of £; recent observations of H/3 hydrogen line emission from interstellar 
H I (Reynold 1971) yield a direct measure of the recombination rate, and the 
corresponding ionization rate in the regions examined could be as high as 10“ 14 
(sec H atom) -1 . 

Since essentially all the gas in the Gum Nebula has been ionized by the 
supernova explosion, the average rate of ionization in interstellar hydrogen by 
such events would be given by 


i <V = f SN N . (18 > 

where <n„ > is the mean interstellar hydrogen density and f is the average 
supernova frequency per unit volume in the galaxy. If supernovae occur at a rate 
of 1 per 100 years in a volume of 4 x 10 66 cm 3 , f ~ 8 x 10" 77 cm" 3 sec " 1 . 

For N - 2 x 10 62 electrons (Brandt et al. 1971), we obtain £<n > ~ 1.6 x 10 -14 

C H 
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sec -1 cm" 3 . Radio observations at 21 cm indicate a mean value for the neutral 
hydrogen density in the galactic plane of about 0.7 cm -3 (Spitzer 1968); however, 
optical depth effects probably raise < n H > to ~1 cm" 3 . A further increase by 
about a factor of 2 in the value of < n H > may be due to the presence of molecu- 
lar hydrogen in dense clouds. Hence the mean galactic value of £ as indicated 
by the parameters deduced for the Vela X supernova is about 8 x 10" 1S sec" 1 
(H atom) -1 . This ionization rate could be reduced if the frequency of supernova 
explosions of the type that produced the Gum Nebula is lower than the assumed 
rate of 1 per 100 years or if the actual free electron content of the Gum Nebula 
is less than that estimated by Brandt et al. (1971). Furthermore, even if the 
average supernova output is similar to that of Vela X but occurs in a region with 
larger magnetic fields, so that the particles propagate over a shorter distance 
before they lose their energy, the ionized volume produced by the supernova 
would be smaller than that of the Gum Nebula. In this case, the total ionization 
per supernova would be reduced since a larger fraction of the total energy would 
go into heating. Also, those supernovae which occur at more than about 100 pc 
away from the galactic plane ( e.g. : the Crab) will produce energetic particles 
which may escape from the disk without producing significant ionization. 

In order for H II regions such as the Gum Nebula to merge into the H I of 
the interstellar medium, the recombination and cooling times of the clouds 
should be much less than the time between supernova explosions t SN in the 
volume of the nebula. For f SN = 8 x 10" 77 cm -3 sec"i and V = 1.35 x 10 63 cm3, 
t sn ~ 3 x 10 s years. 

The recombination time of hydrogen is approximately t rec - 130 T°- 7 /n e 
(Bates and Dalgamo 1962). For n e = 10 cm" 3 and T = 5 x 10 4 °K, t rec -2.5 x 10 4 
years. Since the recombination rate increases with decreasing temperature, 
both t and t „ are much smaller than t cv , so that the ionized clouds in the 
Gum Nebula will rapidly merge into the H I of the interstellar medium. 

In the intercloud medium, for a density < n e > - 0.16 cm -3 , the temperature 
will not decrease appreciably below 10 4 °K in 3 x 10 5 years. The recombination 
time then becomes comparable to t SN so that the degree of ionization is approxi- 
mately constant as a function of time. The fact that the inferred degree of ioni- 
zation of the intercloud medium probably does not exceed about 10% (Hjellming, 
Gordon and Gordon 1969) is consistent with a cosmic- ray model for the ioniza- 
tion of the Gum Nebula, since as discussed above, the energetic charged par- 
ticles from the supernova will deposit most of their energy in clouds and not in 
the intercloud medium. 
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DISCUSSION 


R. McCray: Why did you use a density of 10 cm' 3 for the filaments ? The low 
frequency radio measurements refer to the general mass of the nebula, 
so the temperature reported by Alexander does not apply to the dense 
regions. 

R. Ramaty: Better H-alpha observations are needed. Then we can determine if 
there is an emission measure of at least 600 pc cm -6 over the entire 
beam of the RAE-1 instrument, or whether that brightness (the lowest 
value recorded by Gum) is very localized. 

H. M. Johnson:, It is not clear how Gum determined the emission measure. 

B. J. Bok: I think that the emission measure will probably come out to be of 
order 400 pc cm -6 . 

J. C. Brandt: I want to emphasize that our estimate of the number of electrons 
and the ionization energy does not depend on the emission measure, but 
only on the dispersion measure. 

D. P. Cox: On the particle ionization theory, would the Gum Nebula be gradually 
getting larger ? 

Ramaty: It might; it depends on how many ionizing nuclei are left. If you as- 
sume a simple model with spherical clouds and a dumpiness factor of 
60 or 65, then the distance that one can travel in a straight line before en- 
countering a cloud is equal to the radius of the cloud times the dumpiness 
factor. So there is a good chance that the average nucleus has already 
hit a cloud and lost its energy. 

Cox: But the excess X-ray luminosity requires that the ionizing nuclei are 
still present. 

Ramaty: The excess luminosity may not originate in the Nebula. It may be due 
to discrete sources or to a general galactic background. It might be due 
to charge exchange, if there is still some neutral hydrogen present. A 
search for gamma rays is probably the best test for whether the cosmic 
ray nuclei are still present in the nebula. 


109 



